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Abstract 


The radio light curve of J1415+1320 (PKS 1413+135) shows time-symmetric and recurring U-shaped features 
across the centimeter-wave and millimeter-wave bands. The symmetry of these features points to lensing by an 
intervening object as the cause. U-shaped events in radio light curves in the centimeter-wave band have previously 
been attributed to Extreme scattering events (ESE). ESEs are thought to be the result of lensing by compact plasma 
structures in the Galactic interstellar medium, but the precise nature of these plasma structures remains unknown. 
Since the strength of a plasma lens evolves with wavelength À as X, the presence of correlated variations at over a 
wide wavelength range casts doubt on the canonical ESE interpretation for J1415+1320. In this paper, we critically 
examine the evidence for plasma lensing in J1415+1320. We compute limits on the lensing strength and the 
associated free—free opacity of the putative plasma lenses. We compare the observed and model ESE light curves, 
and also derive a lower limit on the lens distance based on the effects of parallax due to the Earth’s orbit around the 
Sun. We conclude that plasma lensing is not a viable interpretation for J1415+1320’s light curves and that 
symmetric U-shaped features in the radio light curves of extragalactic sources do not present prima facie evidence 
for ESEs. The methodology presented here is generic enough to be applicable to any plasma-lensing candidate. 
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1. Introduction 


The 40 m telescope at the Owens Valley Radio Observatory 
(OVRO) has been monitoring about 1800 blazars at 15 GHz, 
with a cadence of twice per week, since 2008, in support of 
Fermi-GST (Richards et al. 2011). The OVRO sample consists 
of 20.1 Jy flat spectrum radio quasars (FSRQ). The radio 
emission from FSRQs in the centimeter-wave band is expected 
to be dominated by a compact core. The compactness makes 
the OVRO monitoring program an excellent resource to study 
propagation effects such as lensing. 

Lensing is a result of ray deflection caused by spatial 
fluctuations in the refractive index of the medium of 
propagation. In the context of interstellar ray propagation, 
electron density fluctuations in interstellar plasma provide the 
refractive index variation. The natural length scale probed by 
this effect is given by the Fresnel scale, which for Galactic 
plasma refraction is of the order of ~0.1 au. The study of 
plasma refraction is therefore a unique probe of the interstellar 
medium on such small spatial scales. 

Lensing typically leads to multiple (de-)magnified images of 
the lensed source. Even if individual images cannot be 
resolved, the sum total flux-density variations due to the 
relative motion in the source—lens—observer system may be 
observed in light curves. The OVRO survey with its low, 
arcmin-scale angular resolution falls into this category. 
Separating lensing induced variability from intrinsic variability 


galaxies: jets — ISM: structure — scattering 


is difficult. Lensing, however, must, on average, yield time- 
symmetric features. We undertook a manual search for such 
symmetric events in the OVRO archive. Our search revealed 
several lensing candidates, but one source, J1415+1320 (PKS 
1413+135), stood out for several reasons. J1415+1320 shows 
recurring, achromatic, and highly symmetric U-shaped features 
in its light curve, which are unlike anything reported so far (see 
Figure 1 and Section 1.2 below). A collation of multi- 
frequency data from the Metsdhovi radio observatory, the 
Submillimeter Array (SMA), and OVRO shows that the 
U-shaped events persist even in the millimeter-wave band. 


1.1. The Mystery of U-shaped Events 


Time symmetric U-shaped events in the light curves of 
compact extragalactic sources were first reported by Fiedler 
et al. (1987), who named them extreme scattering events 
(ESE). The high degree of symmetry and chromaticity 
(frequency dependence) of such events led Fiedler et al. 
(1987) to interpret them as lensing by intervening plasma 
structures. The plasma lenses were inferred to possess high 
column densities (N, ~ 10!7 cm~ D) over relatively small 
transverse extents (S10 au). Such lenses will be highly over- 
pressurized (density of ne ~ 103 cm™°), and therefore cannot 
exist in pressure balance with the ambient interstellar medium 
(ne ~ 0.03cm~). To mitigate this problem, two classes of 
models were put forth. In the first class of models, the lenses 
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Figure 1. Multi-frequency multi-decade radio light curve of J1415+1320 from various monitoring programs. Flux density offsets have been added to the light curves 


for clarity. The peculiar U-shaped events have been highlighted. 


are ionized sheaths of self-gravitating sub-stellar objects 
(Walker & Wardle 1998). Such objects will, however, have 
to form a considerable fraction of Galactic dark matter if they 
are to generate the observed occultation rates (Walker & 
Wardle 1998). The second class of models invokes highly 
elongated plasma sheets seen edge-on, such that the required 
column density can be obtained at modest over-densities 
(Goldreich & Sridhar 2006; Pen & King 2012). Whether 
plasma sheets with such large axial ratios indeed exist remains 
unknown, and ESEs have remained an unsolved mystery for 
over three decades. 


1.2. The Peculiar Light Curve of J1415+1320 


The incidence of ESEs at 2.7GHz is about 7 x 10-3 
source ! year! (Fiedler et al. 1987). The incidence of 
U-shaped events in the 15GHz OVRO survey is about 
107° source! year | (Vedantham et al. 2017), which is 
roughly an order of magnitude smaller. The large difference in 
the two rates is expected due to the à dependence of the 
plasma refractive index. In any case, the repeating nature of 
U-shaped events in J1415+1320 is unusual. Consider for 
instance, the 2.7 GHz ESE rate. Based on this, J1415+1320 
should undergo an ESE every 140 years if the arrival of ESEs is 
a random Poisson process. We instead observe four U-shaped 
events in about 25 years. This clustering is even more 
significant if we consider the lower incidence of U-shaped 
events from the 15 GHz OVRO survey. 

The light curve of J1415+1320 is unusual even if we set 
aside the above statistical argument. The strength of a plasma 
lens evolves with wavelength as à, and ESEs have so far only 
been seen at Acem Z 2 (Fiedler et al. 1987; Pushkarev et al. 


2013; Bannister et al. 2016). J1415+1320 instead shows 
correlated ESE-like variability down to Acem = 0.1. First, this 
implies a particularly strong and compact plasma lens, further 
compounding the ESE mystery. Next, the observed magnitude 
of (de-)magnification is nearly achromatic over a factor of 20 in 
wavelength, which brings into question the very hypothesis of 
plasma lensing. These inconsistencies prompted us to take a 
fresh look at the viability and consequences of the plasma- 
lensing hypothesis—the principal aim of this paper. 


1.3. Outline of the Paper 


To place the peculiar light curves of J1415+ 1320 in context, 
we describe the properties of the source in Section 2. We then 
discuss the feasibility and implications of plasma lensing in 
Section 3. We end with a discussion of our main findings in 
Section 4. We note that though Section 3 is geared toward 
explaining the peculiar light curve of J1415+1320, most 
aspects of our analysis are generic enough to be applicable to 
any plasma-lensing candidate in the radio band. As such, this 
section is a case study in plasma lensing of extragalactic radio 
sources. 


2. Source Properties and Lensing Geometry 


We now briefly summarize the properties of J1415+1320 
and the plausible lensing geometries. 


2.1. Optical Association 


J1415+1320 is classified as a blazar (on-axis AGN) of the 
BL-Lac type (Beichman et al. 1981). The core of the radio 
source is co-incident with the optical light centroid of an 
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edge-on spiral galaxy at z ~ 0.24 (Carilli et al. 1992, HI 
absorption data), albeit an angular offset of 13 mas (50 pc) has 
been reported based on the isophotal analysis of Hubble images 
(Perlman et al. 2002). In addition, bright radio sources such as 
J1415+1320 are almost exclusively associated with giant 
elliptical galaxies. These factors have led to suggestions that 
the radio source is a background object unrelated to the spiral 
(Stocke et al. 1992; Perlman et al. 1994). If so, then the absence 
of multiple images due to gravitational lensing by the bulge of 
the spiral constrains the radio source to be no further than a 
redshift of z < 0.5 (Lamer et al. 1999), 

Optical spectroscopy of the spiral galaxy by Vedantham 
et al. (2017) shows several lines consistent with a redshift of 
z = 0.247 as expected. There is no evidence for emission 
from another redshift. No broad lines from the AGN 
were detected, most likely due to extinction caused by dust- 
lanes in the spiral. Finally, no excess line emission around 
6563 + 50 A (Galactic Ha), is seen up to a 3o limit of 
2.5 mJy. 


2.2. Parsec-scale Radio Structure 


VLBI imaging of J1415+1320 reveals a relatively young 
radio source with a bright core and a two sided jet-lobe structure 
straddling the core (Perlman et al. 1996). Based on this 
morphology, J1415+1320 belongs to the family of compact 
symmetric objects (CSO). The bulk of the observed radio flux 
shortward of about 2cm comes from the bright unresolved 
core (Perlman et al. 1996). At the observed source flux-density 
of ~ 1 Jy, assuming a jet Lorentz factor of y= 1 and a 
brightness temperature of T, = 5 x 10!° K, we expect the core 
to have an angular size of o,(A) = 90Acm pas. Here the source 
is assumed to be a Gaussian with FWHM of 2.3550,. The core 
remains unresolved in 15 and 43 GHz VLBI images (Perlman 
et al. 2002), which independently yield o,(15 GHz) < 200 uas 
and o,(43 GHz) < 85 yas, respectively. The core, however, 
may internally consist of compact pas scale components whose 
brightness temperatures may be as high as J, ~ 10! K, as has 
been observed in other sources (Johnson et al. 2016). 


2.3. Angular Broadening 


Angular broadening due to scattering in interstellar plasma 
will have an additional influence on the apparent size of the 
source. Since the source is at a high Galactic latitude 
(b ~ 65deg), angular broadening in the Milky Way is 
expected to be very small (a, S 1 pas at 15 GHz; Cordes & 
Lazio 2002). Angular broadening in the interstellar medium 
of the putative spiral host, however, could be significantly 
larger due to its edge-on configuration and because the radio 
source is seen through a giant molecular cloud in the spiral 
(Wiklind & Combes 1997). Very little is known about the 
turbulent properties of interstellar plasma in other galaxies. 
We may, however, obtain rough estimates by assuming that 
the spiral has a “Milky Way like” interstellar medium, and by 
appropriately scaling the value of angular broadening 
expected along an equivalent sight line in the Milky Way. 
To account for the angular offset between the radio source 
and the optical light centroid of the spiral, we consider 
a Galactic longitude and latitude of Odeg and 0.5 deg 
respectively. The expected angular broadening for such a 
Galactic sight line is about 0.5 mas at 15 GHz (Cordes & 
Lazio 2002). In addition, the magnitude of angular 
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broadening must be scaled by the geometric “lever arm” 
as!! Dj; /D;. If the radio source is in the spiral, the lever arm is 
small and of the order of Dis /D; ~ 10 kpc/824 Mpc ~ 1075. 
The expected angular broadening in this geometry is at 
a sub-yuas level and can be neglected. If instead, the radio 
source is behind the spiral, for Dıs = 300 Mpc, the expected 
angular broadening is 0.5 mas x D),/D, ~ 180 yas. This 
value is marginally inconsistent with multi-frequency VLBI 
maps of J141+1320 (Perlman et al. 1996). We note, however, 
that the interstellar medium is known to be “patchy” and 
the line of sight to J1415+1320 may pass through a region 
of significantly lower turbulence. Hence, the above cannot 
be admitted as evidence against the background-source 
hypothesis. 


2.4. Plausible Lensing Geometries 


Given the above known properties of J1415+1320, we 
consider four plausible lensing geometries arranged in increas- 
ing distance between the Earth and the lens (see Figure 2 for a 
cartoon depiction). 


1. The lens is a Galactic object. Since D, >> Dı, whether the 
radio source is associated with the spiral or not is 
inconsequential. The high galactic latitude of J1415 
+1320 (b ~ 65 deg) implies that the lens distance is 
likely not larger than the scale height of the Galactic disk: 
D; < 0.3 kpc, failing which the lenses will have to be an 
exotic halo population. 

2. The lens is in some intervening system that is unrelated to 
both the spiral and the Milky Way. Here the association 
of the radio source with the spiral has only a moderate 
effect on lensing parameters. Since no intervening galaxy 
has been detected, this model essentially requires the lens 
to be associated with an optically faint object such as a 
distant isolated globular cluster or a faint dwarf galaxy. 

3. Both the source and lens are in the spiral galaxy at z = 0.24. 
This gives, D; ~ D, = 824 Mpc, and due to the edge-on 
configuration of the spiral, we assume Dı, < 30 kpc. 

4. The lens resides in the spiral (D; = 824 Mpc), but the 
radio source is an unrelated background object. As argued 
in Section 2.1, z < 0.5, which gives Dı, X 610 Mpc. 


~ 


We now explore the feasibility of, and evidence for, the 
plasma-lensing hypothesis for the above geometries. 


3. Plasma Lensing 


The refractive index of plasma is less than unity. A plasma 
lens therefore diverges light rays. This leads to de-magnifica- 
tion of occulted sources when the observer—lens—source all lie 
on a straight line. At small displacements, curvature in the lens 
profile leads to magnification events that bracket the de- 
magnification trough. This leads to the characteristic U-shaped 
profile seen in extreme scattering events. A cartoon depiction 
of this phenomenon is shown in Figure 3. We now determine 
the minimum lensing strength required to explain the light 
curve of J1415+1320. 


11 : ; 
We use the standard notation wherein D;, D,, and D;, are the observer-lens, 
observer—source, and lens—source angular diameter distances respectively. 
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Figure 2. Cartoon showing the four lensing geometries under consideration and typical distances in the observer—lens—source system. 
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Figure 3. Cartoon depiction (not to scale) of the effect of a diverging plasma 
lens on a background source. The black circle is the Earth, and the gray circle is 
the lens. Relative motion in the source—lens—observer system leads to 
symmetric magnification events that bracket a de-magnification trough as 
shown in the U-shaped curve on the left. 


3.1. Minimum Lensing Strength 


The strength of a lens depends on its intrinsic properties as 
well as the lensing geometry. These two properties can be 
combined into a dimensionless “lensing potential” Y(@), where 
0 is the two-dimensional image position. The source position 8 
corresponding to image position @ is then given by the lensing 
equation: 


B = 0 — Vot (0). (1) 


For a given 6, the lens equation may have more than one 
solution in 0®—a case that corresponds to multiple images of the 


background source. The magnification of each image is given 
by the inverse of the determinant of the Jacobian matrix 
evaluated at the image location: 


1 
~ |det(Ay|’ 


Equations (1) and (2) can be used to determine the location and 
magnification of images for a given source position 6 and 
lensing potential ®(@). See Appendix A for a demonstration. 

The lensing potential for a generic plasma lens is given by 
(Tuntsov et al. 2016) 


ll A = 03/00. (2) 


WA) x Dis 
D:D 


ls 


Tre XN.(0), (3) 


where rẹ is the classical electron radius and N,(0) is the 
electron column density within the lens. Lensing requires two 
conditions to be satisfied: (1) the source angular size should be 
comparable to or less than the lens angular size and (2) the lens 
should possess sufficient strength to deflect light rays by angles 
comparable to the angular size of the lens. Consider for 
instance, a Gaussian lens of the form 


lar 
Y0) = a(Ajexp| ——, |: (4) 
207 


The lens has an effective solid angle of Q; = 2ro? and the 
strength parameter 
1 Ds 


2n07 DD, 


a(A) = Nete X, (5) 
where N, = N,(0) is the transverse gradient of the lensing 
potential on angular scales comparable to the size of the lens 
itself. If a < 1, the effect of the lens is negligible. For a < 1, 
(de-)magnification is moderate, and the lens equation has a 
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Figure 4. Light curves for a point-like source occulted by a Gaussian plasma 
lens of unit size (q = 1). 


single valued solution (no multiple images). a = 1 leads to 
multiple images, and large (de-)magnification. Figure 4 shows 
the light curves for a point-like source occulted by a plasma 
lens with a column density profile given by Equation (4) for 
a = 0.1, 1, and 10. For a > 1, two pairs of caustic are visible. 
We have simulated several other (non-Gaussian) lens profiles 
and find that while the inner pair of caustics are ubiquitous, the 
presence of an outer caustic pair is not a generic property. 

J1415+1320 shows large variations in magnification over a 
factor of 20 in wavelength or a factor of 400 in lensing 
potential. A lower bound on the lensing potential is set by the 
requirement for lensing to have a significant effect on ray 
propagation at the shortest wavelength of Ag = 0.1. More 
formally, we require a(A = 0.1 cm) > 1, or a(A = 2 cm) > 
400. Normalizing to a nominal lens size’? of a = 100 pas, 
Equation (5) yields the constraint 


2 
aa cm? (6) 


Equation (6) places a lower limit on the strength of the plasma 
lens that is conservative due to two reasons. (1) In our analysis 
thus far, we have assumed a point-like source. For sources 
comparable to or larger than the lens angular size o,, and the 
lensing strength, a must be further increased to generate similar 
light curves (see Figure 8 for a demonstration). (2) Larger 
lensing strengths can yield higher order caustic structure that 
result in multiple U-shaped events per lens crossing—a 
scenario that may still be consistent with the data. See, for 
instance, the a = 10 curve in Figure 4, which has one U-shaped 
event on either size of the lens centroid. 

We now explore the consequences of the implied column 
density from Equation (6). We discuss the consistency of the 
observed light curves with a Gaussian lens model constrained 
by Equation (6) in Section 3.4. 


3.2. Free—Free Opacity 
Table 1 summarizes the direct consequence of the limits set 
in Equation (6). We have assumed an isotropic lens—the lens 


12 This will make the lens comparable in size to the upper limits on the source 
at 15 GHz. 
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has the same line-of-sight size as its transverse size, and 
obtained the electron density, ne, which in turn yields the free— 
free optical depth, Te œ f ne di (details in Appendix C). The 
free—free optical depth must be significantly smaller than unity, 
failing which radio photons cannot freely pass through the lens 
and the magnification events from plasma lensing cannot be 
obtained. 

All geometries with an isotropic extragalactic lens may be 
rejected based on their high free—-free opacity. Next, we 
consider anisotropic sheet-like geometries viewed edge-on. For 
such lenses, the free—free opacity from Table 1 must be reduced 
by a factor equal to the axial ratio: ¢ = depth/width of the lens. 
All extragalactic lenses still require unrealistic axial ratios of 
>10? to allow passage of 15 GHz photons. 

The free—free opacity value scales as 0}, and it is therefore 
appealing to rescue extragalactic lens models by significantly 
reducing the lens size to <10 pas. The sources being lensed 
will then have to be compact ~pas scale components in the 
relativistic jet. For instance, lens modeling requires lensing of a 
source component of about 600 mJy in addition to an unlensed 
component (details in Section 3.4). For a relatively high 
brightness temperature of J, = 10! K, the implied component 
size at A = 2 cm is about o ~ 10 pas, which could be lensed 
by a ~10 pas scale lens. A relativistically moving 600 mJy 
source component can cross four 10 pas scale extragalactic 
lenses separated by ~50 pas, to yield four year-long events 
separated by five years each. The distance traversed by such a 
source component is about 200 jas with respect to a 
presumably unlensed stationary core of comparable brightness. 
In addition, lensing must yield multiple images near caustic 
crossings with characteristic separation at A =2cm of 
~ga = 4 mas. Neither the appearance and disappearance of 
such components around caustic crossings, nor bright 
(~600 mJy) relativistic components have been observed in 
multi-epoch VLBI monitoring of J1415+1320 at A= 2cm 
(Perlman et al. 2002; Lister et al. 2016). 

Scenarios where the four U-shaped events are caused by four 
different 600 mJy source components crossing the same lens 
are also problematic since they require a contrived arrangement 
where the source components remain relatively stable during 
the ESE and “cool” rapidly thereafter to maintain the total flux 
density of the source at the ~1 Jy level as observed. Hence, 
apart from free—free opacity constraints, extragalactic lens 
models have an additional problem in explaining the timing, 
duration, and flux densities of the four ESE-like events. We 
therefore find Galactic lenses at D; = 10 kpc (value justified in 
Section 3.3) to be the most viable ESE scenario among the four 
geometries. The lens crossing in such geometries is determined 
by the relative motion between the Earth and the lens, which 
we explore next. 


3.3. Retrograde Motion of Galactic Lenses 


The symmetry of U-shaped ESEs is a result of rectilinear 
relative motion between the source and an axially symmetric 
lens. The orbit of the Earth around the Sun can perturb an 
otherwise rectilinear trajectory leading to asymmetries in the 
observed light curves. The ESE-like events presented here last 
about 400 days (interval between magnification events), which 
is significantly longer than previously reported ESEs. If the 
lenses are Galactic in origin (as implied by Section 3.2), this 
allows us to place unique constraints on the size of the lens as 
follows. Consider a plasma lens at 1 kpc from the Earth. 
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Table 1 
Table Showing Lower Bounds on Plasma Density and Free—Free Opacity for the Putative Plasma Lenses 
Model Dis D; Ne Lens Size Free—Free Optical Depth (15 GHz) 
Galactic lens 800 Mpc 10 kpc 10'° 5 x 10° lau 3 x 10-3 
Lens in intervening dwarf 700 Mpc 400 Mpc 10” 9 x 10° 0.2 pe 5 x 10° 
Lens and source in putative host 30 kpc 800 Mpc 10” 1.1 x 10'° 0.4 pe 1.7 x 10"! 
Background source; lens in putative host 400 Mpc 800 Mpc 10” 1.4 x 10° 0.5 pe 2x 10° 


Note. Rows correspond to different lensing geometries (see Figure 2). A lens size of a = 100 pas has been assumed. Free-free opacity is proportional to ø}. 
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Figure 5. Effect of the annual parallax on lensing light curves for a Galactic plasma lens. Left panel shows the data (2015 event) and model light curves for varying 
magnitudes of this parallactic effect. Right panel shows the likelihood function (1, 2, and 3ø contours over-plotted) for varying magnitude of the parallactic effect. The 


x and y axes lie along the ecliptic plane with x-axis along zero ecliptic longitude. 


Because the parallax due to Earth’s annual motion around the 
Sun is 1 mas at 1 kpc, a lens that is smaller than 1 mas must 
yield a light curve with significant asymmetry. The magnitude 
of this effect is inversely proportional to the product of the lens 
angular size and distance, oD. To test for the presence of this 
“annual parallax” in the data, we simulated light curves for 
various values of oD; and evaluated the residues between the 
data and our model within the U-shaped trough at 15 GHz. We 
assigned the lowest least square error value thus obtained as the 
“noise” in our data and computed the likelihood of all oD, 
choices by assuming Gaussian noise statistics. Anticipating a 
slight degeneracy between oD, and the exact epoch of 
occultation, we marginalized the likelihood over the latter. 
The resulting likelihood function along with the 1, 2, and 3a 
contours is shown in Figure 5 (right panel). 

Using Figure 5, the product of lens distance and lens size 
may be constrained to be œD; > 2.0 mas kpc. D; is unknown, 
but the lensing scenario requires œ ~ o; at 15 GHz (see 
Section 3.4). The best observational constraint yields 
Os < 200 pas (see Section 2.2). The corresponding constraint 
on the lens distance, therefore, is D; > 10 kpc. Since the source 
is at a high Galactic latitude (b ~ 65 deg), these values are 
clearly incompatible with the plasma lens being in the Galactic 
disk where most of the interstellar plasma lies. 


3.4. Inconsistency of Light Curves with Plasma Lensing 


Fitting of model ESE light curves to data is complicated by 
the fact that almost identical light curves can be obtained by 
increasing the lensing strength, while appropriately increasing 
the source size (see Figure 8 for a demonstration). We therefore 
proceed with the most conservative assumption: the source is 
point-like at the highest frequency, 234 GHz. We thus obtain 
the smallest value of a(234 GHz) that can account for the 


range of variability seen in the data. This value of a can then be 
scaled to lower frequencies, and the source size adjusted to 
obtain reasonable fits to the data. We found that matching the 
models to the data required (1) a source size that scales with 
wavelength as X (o/a = 1 is required at 15 GHz) and (2) a 
600 mJy component of the source was lensed at 37 and 
15 GHz. A significantly smaller flux density will not yield the 
required depth in the U-shaped trough. If, on the other hand, 
the entire source is lensed, the large de-magnification at 
15 GHz should have made the source virtually undetectable at 
15 GHz for about a year. The resulting model light curves are 
plotted with the data at 15, 37, and 234 GHz in Figure 6 for the 
2010 and 2015 events (upper and lower rows respectively). 

The model provides poor fits to the data outside the de- 
magnification trough. Part of the reason for this may be our pre- 
supposed (and hence arguably inaccurate) lens profile and 
intrinsic source variability. However, the main reason is the 
generic incompatibility of the data with a de-magnification 
trough that is a hallmark of plasma lensing. Consider, for 
instance, the 15 GHz flux densities of the source (see Figure 1) 
at years 2009.5 and 2011.5. A de-magnification should have 
yielded significantly lower flux in 2009.5 as compared to 
2011.5. The data instead show the source to be of comparable 
brightness at the two epochs. 


4. Discussion and Summary 


J1415+1320 shows recurring and symmetric U-shaped 
features in its three decade long light curves at frequencies 
ranging from the centimeter-wave band to the millimeter-wave 
band. The canonical interpretation for such U-shaped events in 
the radio source is lensing by compact interstellar plasma 
structures. We have critically examined the plasma-lensing 
interpretation and arrived at the following conclusions. 
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Figure 6. Data overlaid with model plasma-lensing light curves (Gaussian lens) for the 2010 and 2015 U-shaped events (top and bottom rows respectively). The 
column are for three frequency bands: 15, 37, and 234 GHz spanning a factor of 400 in lensing strength, a. The ratio between source size and lens size, o/os is given 


in the panel labels. 


1. All extragalactic plasma lenses require unrealistically 
high axial ratios (width/depth) of 210° to yield 
acceptable amounts of free—free opacity (Table 1). 

2. The absence of significant asymmetry expected due to the 
Earth’s orbit around the Sun yields a lens distance that is 
well beyond the Galactic disk. This places the putative 
Galactic lenses in a region, which is expected to be 
largely devoid of dense interstellar material (Figure 5). 

3. The U-shaped events seen in J1415+1320 are not well 
described as de-magnification events. The flux-density of 
the source at the bottom of the U-shaped trough is 
comparable to that away from the U-shaped events. For 
this reason, model light curves provide poor fits to the 
data even in the immediate vicinity of the trough 
(Figure 6). 


The above factors compel us to reject the ESE hypothesis for 
J1415+1320. We therefore conclude that symmetric U-shaped 


radio light curves do not present prima facie evidence for 
plasma lensing. This raises two pertinent questions. 

What caused the U-shaped events in J1415+1320? The high 
degree of symmetric of the U-shaped events in J1415+1320 
points to lensing as the cause. The only other lensing 
mechanism known is gravitational lensing, which is achro- 
matic. The viability and consequences of the gravitational 
lensing hypothesis are discussed in a separate paper (Vedan- 
tham et al. 2017). 

How can one ascertain and constrain plasma lensing? The 
present study of J1415+1320 has shown that largely achro- 
matic U-shaped events may be obtained despite plasma 
refraction being a chromatic phenomenon. This is because an 
increase in source size with wavelength can partly offset the 
increase in lensing strength with wavelength to yield a 
U-shaped event with comparable levels of de-magnification 
over a wide wavelength range (Figure 8). It is the incompat- 
ibility of data with model ESE light curves over a longer term, 
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Figure 7. Image magnification and locations for a Gaussian plasma lens with strength parameter œa = 10. Top and bottom left panels show magnification of image 
components and the sum total magnification respectively. The right panel shows the location of the image components. Black curves correspond to source position that 
yields a single image, and red, green, and blue curves correspond to cases with three images. 


among other constraints, that enables us to rule out the ESE 
hypothesis for J1415+1320. In the absence of longer-term 
monitoring data, the cessation of ESE-like variability at some 
high frequency may be taken as evidence in favor of the ESE 
hypothesis. 
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Appendix A 
Image Location and Magnification 


Figure 7 shows the magnification and location of images 
formed when a point-like source crosses a Gaussian lens 
described by Equation (4), with strength parameter a = 10. For 
simplicity, we have considered a one-dimensional lens. The 


right panel plots the roots of the lens equation, 0; for varying 
source position 3. The lens equation has either one root (i = 0) 
or three roots (i € [0, 1, 2]) depending on source position. The 
magnification of each image (top left panel) is given by 
u; = (03/00) '|g~o,. The light curves presented in this paper 
do not resolve the individual images and are therefore only 
sensitive to the sum total magnification (bottom left panel). The 
central de-magnification trough around 8 = 0 is bracketed by 
high magnification caustics. This yields the characteristic 
U-shaped ESE light curve. Successive U-shaped features in 
light curves may be caused by additional sources crossing the 
same lens or a single source crossing multiple lenses in 
succession. 


Appendix B 
Lens Strength—-Source Size Ambiguity 


Many reported ESEs have either been observed at a single 
frequency channel or totally uncorrelated variations are seen 
between two frequency bands. This leads to considerable 
problems in interpreting the light curves. Several parameters 
such as lens strength, source size relative to the lens, and the 
flux densities of the lensed and unlensed components must be 
fit to observed light curves, leading to degeneracies between 
parameters. In particular, we find a strong degeneracy between 
the source size relative to the lens and the lensing strength. As 
an illustration, consider the well sampled 15 GHz ESE during 
2009. In Figure 8, we have shown two model light curves 
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Figure 8. a—o, degeneracy (a = lensing strength; o, = source size in units of lens size). The two models shown in the two panels both describe the data points equally 
well. It is therefore not possible to infer lensing parameters at a single frequency, or a closely spaced sub-set of frequencies, without explicit knowledge of the 


source size. 


whose lensing strengths a, vary by a factor >500 and yet, the 
curves provide a reasonable fit to the data points. Figure 8 
shows that one cannot determine the lensing strength with light 
curves at a single frequency. Rather, only a lower limit on the 
lensing strength can be placed (assuming the source is 
significantly smaller than the lens). In our case, if only the 
15 GHz light curves were available, we would only conclude 
that a = 0.7 


Appendix C 
Free—Free Opacity 


Let the transverse size of the lens be 6), distance to the lens be 
D,, the electron column density be N,, and let the lens have a 
line-of-sight extent p times that of its transverse extent. The 
electron volume density is then ne = N,/(D,9,C). The emission 
measure (line-of-sight integral of n2) is then EM = Ne / (D10). 

The free-free opacity at radio frequency is approximately 
given by 


-2.1 
Te = 3.28 x 1077 (E) ED ; (7) 
GHz cm~° pe 


where T = 1047; is the gas temperature in kelvin, and we 
assume 7, = 1. 
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